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damping mechanism whose experimental observability could be relevant to understand decoherence 
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I. INTRODUCTION 



> 

Since Schrodinger, both theoretical and experimental efforts have been devoted to understand whether quantum 
mechanical laws still work at the macroscopic level 0. In particular, many avenues are currently investigated to 
£f) create coherent superpositions of macroscopically distinguishable states starting either from microscopic degrees of 
<0 freedom, as in the case of photons electrons ||, neutrons @, atoms j?],|fl and molecules ||, or macroscopic 

dnes, for instance by using superconducting circuits j[o|. The common goal of these attempts is to grasp how 
^^the coherence properties, the crucial role of which is well established in the microworld, are washed out in the 
transition to the macroscopic realm, through what is commonly termed decoherence |11 This last concept, 
Ojnvoked to explain why superpositions of macroscopic states are so fragile to make their actual observation far from 
+L> being a trivial task, is present whenever the system under consideration is coupled to the outside world, turning its 
Ch evolution into an irreversible one. Among the possible sources of decoherence, the one attributable to the act of the 
^measurement deserves special attention since, being any physical inquiry on the system performed via the coupling 
^4° an external probe, it is in principle unavoidable. In previous papers the measured induced decoherence has been 
• analyzed in the case of microscopic systems, namely atoms confined either in optogravitational cavities | fl3| or in 
^ electromagnetic traps Jl4[ . Observable quantities turn out to be modified giving rise to peculiar transient phenomena. 

In this work we extend the analysis to the dynamics of macroscopic quantum systems such as the Josephson devices 
^ which, as emphasized by Leggett |L0]], are promising candidates to study the coherence-decoherence transition by 
Ci virtue of both their low operating temperature and the small intrinsic dissipations. In Section II we introduce the 
framework of measurement quantum mechanics and apply it to the measurement of energy in a radiofrequency 
superconducting quantum interference device (rf-SQUID). Decoherence signatures arising from the measurement 
process in the dynamics of the average magnetic flux are discussed in Section III. Final remarks on the feasibility and 
the differences from the microscopic systems complete the paper. 



II. MEASUREMENT QUANTUM MECHANICS OF A SQUID RING 



The system that will be considered hereafter is a superconducting closed loop interrupted by a weak link of the 
Josephson type, in short an rf-SQUID. Neglec ting dissipation mechanisms, its quantum mechanical description can 
be modeled through the effective Hamiltonian |15[ 
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Pi „ h 2 d 2 - <b ext ) 2 I $n / $ \ 

where C is the weak link capacitance, L the ring inductance, I c the critical current of the junction, & ext is the external 
magnetic flux threading the ring and the flux quantum <£>o = Ti/2e » 2.07 • ICU 15 Wb. In the circuit description ([!]), 
the total magnetic flux <t and the total displacement flux (with units of electric charge) threading across the weak link 
capacitor, P$ = —iTid/d<5>, play the role of generalized coordinate and momentum respectively. Notice that C acts as 
the effective mass of the whole ring, treated as a single macroscopic quantum object evolving in $-space under the 
potential V($>). If the adimensional variable x = ($ — $ ea:t )/$o is introduced and a fixed value & ext of the external 
flux is chosen fulfilling the relation Q ext /&o = (n+1/2), with n integer to ensure ^-symmetry of the effective potential 
around $ = $ ea:t , the latter can be rewritten as 

y (x) = |L 2 + ^cos(27nr). (2) 

ZL ITT 

According to (^), monostable or multistable dynamics arises depending on the relative amplitude of the parabolic 
and periodic terms. In particular, a bistable regime is obtained if 1 < p < 5ir/2, with p — 2ir LI c /§o. Under such 
conditions, the double well potential (||) can be approximated by a polinomial expansion, 

V(x) * V - \x 2 + ^ , M ,A>0, (3) 

being the parameters fj,, A, Vb related to the physical ones via the relationships: 

/i = 27t/ c $ - $l/L , 

A =4^ 3 / c $ /3, (4) 
Vb = /c$o/2tt = 3A/8tt 4 . 

The quartic potential approximation (|3|), possessing two local minima at x± = iy^u/A separated by a barrier of 
height AU = /i 2 /4A, has been widely exploited to model the macroscopic quantum tunneling of the flux through the 
junction jl(|. One can easily check by inverting (^) that only the range of values /i < 3A/27T 2 is allowed in order the 
bistable potential (^) corresponds to a physical SQUID. 

Let us now include the effect of a measurement process in the description of the SQUID ring. We will focus the 
attention on a situation where the measurement most clearly possesses a dynamical content, namely the continuous 
monitoring of a generic system observable A. Significant progresses in modeling such a situation have been made 
in recent years through different approaches, all of which basically recognize that a measured system is not isolated 
but in interaction with a measurement apparatus. A detailed account of these approaches, together with the unified 
picture in which all of them can be incorporated, has been given elsewhere |T^]. Few remarks are in order to apply this 
general scheme to our particular problem at hand. Due to the fact that SQUIDs are macroscopic condensed matter 
objects and therefore two indiscernible rings (let alone an ensemble) cannot be constructed, measurement procedures 
are constrained to deal with just one SQUID and physical information emerges from averaging individual histories 
with the system each time prepared in the same initial macroscopic state. In the language of quantum measurement 
theory, such a procedure translates in the so-called selective a priori measurement scheme, whereby the instantaneous 
result of the measured observable corresponds to a particular, though a priori unknown, selection of the possible 
outcomes. Accordingly, the deterministic Schrodinger evolution for the state vector is replaced by a stochastic one, 
where single realizations closely mimic the fluctuating output of each experimental run. Two equivalent approaches 
may be used for unravelling the master equation, either a Quantum State Diffusion picture where stochasticity acts 
continuously |p"7j| , or the so-called Quantum Jumps approach resulting when stochasticity is chosen to act occasionally 
fliSfl . Nonselective measurement results are recovered by averaging the selective outcomes according to their probability 
distribution. On the other hand, a different description of the dynamics in the presence of the measurement, although 
not viable in practice when no ensemble is at disposal, has been proven to be equivalent to the previous one insofar 
as averaged (hence nonselective) predictions are required. Being in this case technically much more simpler, it will be 
adopted for the SQUID problem. 

The starting point for a direct description of nonselective measurements is the master equation, derived in the 
framework of open quantum systems, which rules the evolution of the density operator p(t) of the system under 
observation Jl9[ , 

^ = -i[H, m ]-~[A,[A,p(t)]], (5) 
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where H, given in (0), describes the dynamics of the unmeasured closed system (k = 0) and the measured observable 
A is a so-called Lindblad operator representing the influence of the environment on the system. The parameter 
k expresses the strenght of the coupling between the measured system and the meter, that will be assumed to be 
time-independent throughout the derivation fl6|| . Let us restrict the discussion to the case A = H, corresponding 
to a continuous measurement of energy, with k = he- Other situations in which the measurement process has 
been taken into account for a bistable potential have been previously discussed in Q an( i H^l f° r continuous and 
impulsive monitoring respectively, but choosing as measured observable the (possibly generalized) coordinate itself, 
i.e. A = $ in (^). Our choice, which reflects a viewpoint already discussed in fl3| , is instead guided by the demand to 
describe the effect of a measurement of one observable on the dynamics of another one which does not commute with 
the former. This also corresponds to the original philosophy of the quantum Zeno effect as introduced by Misra and 
Sudarshan, whereby measurements of the particle track were expected to affect the lifetime p^| . Moreover, the energy 
observable plays a peculiar role for all the models in which a fundamental diagonalization in the energy eigenstates 
(the so-called intrinsic decoherence) is predicted to occur and for which a formally identical description applies |^3|,Q . 
The nonselective master equation (|^) can be conveniently solved by introducing the representation of the measured 
observable eigenstates, H\n) = E n \n), with n — 0, 1, . . . labelling the quantized energy levels of the Hamiltonian (|l|) 
and density matrix elements p nm (t) = (n\p(t)\m) . Eq. (|5|) is thus rewritten as 



whose general solution is 



,(t) - ~\{En - E m ) Pnm (t) -^-(E n - E m ) 2 Pnm {t) , (6) 



Pnm (t) = exp <j --(£„ - E m )t - — {E n - E m yt \ p nm {0) . (7) 



This allows us to write down the exact density operator evolution for the system undergoing the measurement process 
in the form 

#*) = H/W(0)exp j~J(£ n -E m )t- ^{E n -E m ft\ \n)(m\ , (8) 

nm, ^ ' 

where the quantum decoherence effect accompanying the measurement is manifested as the vanishing of the off- 



diagonal density matrix elements at a rate proportional to the associated energy spacing |12j . According to Eq. (|S|) , 
the evolution of the system has been altered with respect to the closed one, ke = 0; it is important to understand how 
such modifications can manifest in observable properties of the system. Our approach, based on a direct observation 
of the monitored quantum system with energy meter variables disregarded, can be considered as complementary to 
the adiabatic monitoring scheme already proposed in |p5| , where the quantum variables are instead traced out to 
leave a modified reduced evolution of the apparatus. 



III. QUANTUM DAMPING OF THE AVERAGE MAGNETIC FLUX 

As a SQUID observable showing nontrivial dynamical behavior, let us consider the total average magnetic flux 
threading across the ring as a function of time, (<£(£)) = Tr($/5(t)) = ^2 nm p n m(t)(m\<&\n) . Owing to Eq. (g), we get 

=^/w(0)exp i~^(E n - E m )t - ^(E n ~ E m ) 2 t\ (m\^\n) . (9) 

nm ^ ^ 

At variance with the purely oscillatory behavior occurring when no measurement is performed, an exponential damping 
with state-dependent time constants r nm = 2/[nE{E n — E m ) 2 ] arises, for which no classical counterpart can be 
envisaged. After a transient whose duration depends, for a fixed ke, on the Bohr frequencies uj nm = (E n — E m )/h 
picked up by the initial configuration, the average flux dynamics becomes completely inhibited leading, irrespective 
of the initial state, to an average localization around the asymptotic value 

lim ($(*)) =Vp„ n (0)(n|4|n) = ! k e > , (10) 

£— > + OC ' 

n 

where odd parity of the flux operator has been exploited. 
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In order to visualize the phenomenological features arising in the presence of quantum damping, we specialize the 
general description of Eq. (g) by choosing some representative initial SQUID configurations of the system. A picture 
of the double well potential (||) and its lowest energy eigenfunctions is reported in Fig. 1. Firstly, let the SQUID 
be initially prepared in a combination of the two lowest states, \^,t = 0) = l/v2[|0) + = \L), where |0) and 
|1) are the symmetric ground state and the antisymmetric first excited state respectively and the notation \L) has 
been introduced to stress that the dominant concentration of the state is in the left well. We will also consider 
the antisymmetric right-localized superposition state \R), defined by 1 / a/2 [|0) — = \R). Due to their localized 
nature, the states \L), \R) can be differentiated on a macroscopic level by the net left or right origin of the magnetic 
flux (also implying an opposite circulation sense for the superconducting current) and represent therefore a simple 
realization of two macroscopically distinguishable states. By evaluating the initial density matrix p(0) = \L)(L\, a 
damped oscillation is found from Eq. ([)]) for the average flux motion, 

<$(*)> = exp (- ^(M 2 i)^{(0|!.|l)e— «} , (11) 

where parity conservation has been again taken into account and uio = to = (E\ — Eo)/h denotes the angular tunneling 
frequency, related to the height barrier and the zero-point energy Hujq in each well via a WKB estimate [[15[ : 



AU f AC/1 

exn < > 



wo\/t— exp<^ - — \ . (12) 
V nujQ ( Tiluq J 

Transient phenomena in the average flux dynamics are ruled in this case by a time scale of the order 2/[ke(Tiu) 2 ]- 
Some interesting insights on the action of the measurement process can also be gained by depicting the evolution of the 
|L)-state directly in the Li?- representation, which is the natural one for tunneling analysis. This may be accomplished 
by rewriting the relevant projectors in (@) according to the inverse transformation, e.g. |0)(0| = [|L)(L| + \R)(R\ + 
\L)(R\ + \R)(L\]/2, |0)(1| = [\L)(L\ - \R)(R\ - \L)(R\ + \R)(L\}/2 and so on. As a result, the density operator of the 
measured system in the macroscopic Li?-representation can be written as follows: 

„. . ^ / 1 + c- KE(n ^ 2t / 2 cosujt -ie- RE( - n ^ 2t/2 smujt \ , , 

P\ > - y ie - K , E (hu J ) 2 t/2 sinujt 1 _ e -KE(huj) 2 t/2 coaujt J ■ \ > 

The off-diagonal density matrix elements make evident the quantum interference between the macroscopically distin- 
guishable states \L) and \R), implying the existence of a Schrodinger cat state 0. According to (|l3|), the decay of this 
macroscopic coherence represents the off-diagonal effect of the coupling to the meter. In addition, the measurement 
influence extends to the diagonal matrix elements, expressing the probability for left and right localization of the state 
respectively. Asymptotically, the tunneling mechanism becomes frozen, a behavior showing close similarities with the 
result quoted in |2q| , where a dynamical suppression of tunneling in a symmetric double well potential arises as a 
consequence of an external perturbation - a coherent laser field driving transitions to an electronically excited bistable 
state. Of course, the same evolution given in ([□]) has to be recovered using the density representation (|l3|), 

($( t) ) = M^(l + c --(MV2 cosw i) + i^l^^e-^tM 2 ^^^) 

= (0|*|0) + (0|4|l)e- KE(?lw)2t/2 cosc^ . (14) 

A more realistic preparation of initially localized states usually takes place by considering Gaussian flux wavepackets 
of the form 

M * ] = KF exp I J = Wo J^W exp I — 2^- J • (15) 

with the parameters <i> m = x m &o + Q ext and cr$ = $oc K properly chosen to attain both the desired flux localization 
and a suitable initial average energy ((ipo\H \?po) < At/ in the tunneling regime). Starting from 



the nonselective evolution of the average flux results from Eq. (|9|) . By referring to |L3|] for general considerations on the 
Gaussian case still applying here, a quantitative description is made available by resorting to numerical calculations. 
The required energy eigenfunctions and eigenvalues have been computed through a selective relaxation algorithm 
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extensively described in j27J] . A representative evolution arising from a Gaussian state maximally contributed by the 
first four energy eigenstates is depicted in the sequence of Fig. 2, where the effect of an increasing coupling to the 
meter proceeds from 2(a) to 2(f). A more resolved picture for the first three cases follows in Fig. 3. The presence 
of the second doublet of almost degenerate eigenstates has been obtained by centering the Gaussian state around a 
value of magnetic flux intermediate between the ones corresponding to the left maxima of the eigenfunctions (Fig. 1). 
The transient strongly depends on the initial state and exhibits, as already remarked, many decay timescales. When 
ke > 0, after a time r only the oscillations involving levels with energies \E n — E m \ < ^/2/(k,et) will survive. More 
in general, critical values of the measurement coupling constant can be introduced for each oscillation pattern as[| 

Tnm(E n E m ^ h{E n E 7n ^ 

where T nrn = 2~k /uj nm are the periods of the Bohr oscillations. Each critical value rules the nm-subdynamics of the 
SQUID magnetic flux, giving rise to deviations from the oscillatory closed system evolution for ke > K-Zm- Since 
the smaller energy splitting is due to the ground state {\E\ — Eq\ > \E n — E m \ Vn, m), the overall dynamics will be 
affected if ke > k^q 11 '■ By varying ke it is thus possible to observe qualitatively different initial transients, selecting 
for instance the oscillations at the frequencies LU32 and u>io (Figs. 2(b) and 3(b)), or the latter only (Fig. 2(c-f)). 
Further increase in causes the oscillations to disappear at all, Fig. 2(c-f). 

Few remarks are worthwhile. Similar considerations can be in principle repeated for the evolution of any ensem- 
ble average quantity, the explicit solution of the appropriate stochastic Schrodinger equation (or Quantum Jumps 
equation) being required to calculate individual SQUID behavior. Note, however, that the transient nature of the 
expected phenomena would not allow us to exploit time averages of the monitored quantities sampled over a single 
stochastic trajectory, the former being only viable if the observed system is at thermal equilibrium as discussed in [|28| . 
Indeed, having thermal fluctuations been neglected in our case, the same argument given in ]Tij ] can be applied to 
show that, for a SQUID density matrix evolving according to Eq. convergence of individual runs to energy eigen- 
states is expected. Conceptually, the damping mechanism operating in the nonselective ensemble evolution of average 
properties can then be also interpreted as due to the random dephasing between different histories associated to the 
selective evolutions of the same properties - both dephasing and damping being signatures of irreversible dynamics. 
In addition, in order to make the predicted behavior qualitatively similar to the one occurring in practice, the physical 
parameters have to be chosen in such a way that both bistable and quantum tunneling regimes are guaranteed by the 
potential energ y (3 ) . According to the values of /i, A corresponding to the Fig. 1, an ultra small capacitance SQUID 
ring as used in [P9| has been modeled, C ~ I CU 16 F, leading to L ~ I CP 11 H, I c ~ 20 fj,A. For such a configuration, 
which is not too distant from routine design, thermal activation can be made negligible at relatively high operating 
temperatures, e.g. k B T/AU ~ 6 • ICT 3 at T = 4 K. 

IV. DISCUSSION AND CONCLUSIONS 

A quantum damping due to the interaction between a macroscopic quantum system and an informational environ- 
ment has been discussed in the case of a SQUID ring undergoing continuous nonselective measurements of energy. 
This damping has no classical analog and is strongly imprinted by the decoherence process originating in the act of 
measurement. 

The predicted effect can be made observable for instance by measuring the dynamics of the average magnetic flux 
in the SQUID or any other observable which does not commute with the Hamiltonian (Q). Since the system under 
observation is a single macroscopic object the procedure must consist in repeating the cycle of measurements, each 
characterized by the following steps. Firstly, the system is prepared in a flux state, say a Gaussian, centered around 
one of the two minima for instance by properly adjusting the external magnetic flux. At the same time the energy 
meter is turned on. The total magnetic flux is measured after a given duration t by means of an instantaneous von 
Neumann measurement. The cycle is repeated again at constant t allowing to build the magnetic flux distribution 
and therefore the average magnetic flux. The whole procedure is repeated again by sweeping the time t, allowing a 
comparison with the predicted behaviour (^J). Concerning continuous measurements of energy, to our knowledge no 
experimental technique is by now available. The closest achievement has been reported by the Sussex group regarding 
the spectroscopy of quantum mechanical SQUID rings p9[ . The kind of measurement required here is however quite 
different, since a dynamical monitoring of energy is demanded. In other words, even if the knowledge of the energy 



1 For analogous definitions of critical measurement couplings in atomic systems see for instance pAjjL 
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levels is a prerequisite, a monitoring of the populations in each energy eigenstate is also necessary. For states formed 
by the lowest two levels, this could be achieved by coupling the two-level dynamics to a third level in a way similar 
to the one exploited to monitor the electronic populations in the study of atomic quantum Zeno effect [ pQ ]. In the 
SQUID problem, the lack of available schemes does not allow to make quantitative comparison between the predicted 
time decay due to quantum damping and the decays times due to classical sources of decoherence associated to 
intrinsic dissipations |3l) . In addition, if fundamental energy localization is proven to occur in nature, the detection 
of measurement induced energy localization will only be possible if the measurement coupling ke is sufficiently large 
compared to the minimum step time associated to intrinsic decoherence |23| , |24|| . In any case, we guess that its actual 
observability is obtainable for values of the parameters not far from those required to make intrinsically quantum 
phenomena, such as macroscopic quantum coherence J3^] or temporal Bell inequalities [jl0|,^3), detectable. It is also 
worth to point out that the predicted damping can be also implemented in a more quantitative way by considering 
tunneling phenomena at the microscopic level, for instance in the spectroscopy of molecular systems as considered in 
p6[ . On the other hand, in this case the possibility to test decoherence of a macroscopic degree of freedom is lost. 

The relevance of quantum damping in post-modern quantum mechanics is twofold. Firstly, it is an inequivocable 
test of the validity of current models on the quantum-to-classical transition which identify decoherence induced by the 
opening of the quantum system to the external world as the key concept. The possibility to prepare states in such a 
way to produce controllable decay patterns can support evidence for or against quantitative predictions of decoherence. 
Moreover, the decoherence issue in SQUID rings is crucial to exploit superconducting circuits as quantum computers 
. The model presented here can be considered a highly idealized one describing the ultimate source of decoherence 
due to the informational input-output operations performable on superconducting quantum gates. 
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FIG. 1. Plot of the double well potential V(x) - V of Eq. @ for the parameter choice li = 1.80487(1) eV, A = 14.73360(1) 
eV, resulting in local minima at \x±\ = 0.35 and height barrier AU = 0.055274 eV. The four lowest energy normalized 
eigenfunctions <fi n (x) (n = 0, . . . , 3) are also shown together with their energies Eo = —0.0440591 eV, E\ = —0.0440585 eV, 
E 2 = -0.0231600 eV, E a = -0.0230991 eV. 

FIG. 2. Predicted time dependence, in units of the tunneling period Tbi, of the average magnetic flux for different values 
of the measurement coupling relative to the critical values «io lt and Ky^ 1 and for an initial Gaussian wavepacket having 
x m — —0.27 and o x — 0.06. The state is contributed up to w 97% by the two lowest energy doublets, (a) ke = 0; (b) 

in-2^cri£ in - 4,^cri£. / \ in-l^. cr ^ in-3,^crii. / t\ ,„ inherit in - ^ «crit, / \ in2 ,crit „crit 

KB — 10 K, 32 — 10 K W ; (C) HE = 10 K 32 = 10 K 10 | (d) K.E = WK, 32 = 10 K 10 ; (e) HE = 10 ft 32 = K 10 \ 

(f) he = 10 3 k§"* = IOkio*'. The dark band in (a)-(c) around the tunneling motion with luio is due to the high frequency 
secondary oscillations picked up by the initial state. 

FIG. 3. Same as Fig. 2 on a time scale expanded by an order of magnitude for the previous cases (a), (b) and (c). In (a) 
all the Bohr frequencies are present; in (b) only the slower oscillations at lu\o and U032 survive, with UJ32 /low ~ 10~ 2 . The 
component disappears within a tunneling period for a larger value of ke, as in (c). 
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